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ABSTRACT 

Recent observations suggest that Type Iln supernovae (SNe Iln) may exhibit late-time (>100 days) 
infrared (IR) emission from warm dust more than other types of core-collapse SNe. Mid-IR observa- 
tions, which span the peak of the thermal spectral energy distribution, provide useful constraints on 
the properties of the dust and, ultimately, the circumstellar environment, explosion mechanism, and 
progenitor system. Due to the low SN Iln rate (<10% of all core-collapse SNe), few IR observations 
exist for this subclass. The handful of isolated studies, however, show latc-timc IR emission from warm 
dust that, in some cases, extends for five or six years post-discovery. While previous Spitzer /IRAC 
surveys have searched for dust in SNe, none have targeted the Type Iln subclass. This article presents 
results from a warm Spitzer /IRAC survey of the positions of all 68 known SNe Iln within a distance of 
250 Mpc between 1999 and 2008 that have remained unobserved by Spitzer more than 100 days post- 
discovery. The detection of late-time emission from ten targets (~15%) nearly doubles the database 
of existing mid-IR observations of SNe Iln. Although optical spectra show evidence for new dust 
formation in some cases, the data show that in most cases the likely origin of the mid-IR emission is 
pre-existing dust, which is continuously heated by optical emission generated by ongoing circumstellar 
interaction between the forward shock and circumstellar medium. Furthermore, an emerging trend 
suggests that these SNe decline at ^1000-2000 days post-discovery once the forward shock overruns 
the dust shell. The mass-loss rates associated with these dust shells are consistent with luminous blue 
variable (LBV) progenitors. 

Subject headings: circumstellar matter — supernovae: general — supernovae: individual — dust, 
extinction — infrared: stars 



1. INTRODUCTION 

Type Iln supernovae (SNe; see iFilippenkd IT9971 for a 
review) have gained considerable attention over the past 
decade. Representing less than 10% of all core-collaps e 
events (ISmarttl [200a ISmith et alllMlt l!Tet~aII [2^11 . 
this subclass comprises more than 50% of the rare sample 
of SNe observed to exhibit a late-time (>100 days) in- 
frared (IR) excess with respect to the optical blackbody 
counterpa rt, typically suggesting the presence of w arm 
dust (e.g.. lGerardv et al.ll2002l: iFox et al.|[2009ll2010l and 
see Table[T|). Named for their relatively "narrow" (< 700 
km s _1 ) emission lines arising from a dense, slo wly mov- 
ing, p re-existing circumstellar medium (CSM) (|Schlegell 
I1990T ) . these SNe Iln may be more prone to have warm 
dust than any other SN type since (1) dust is able to form 
in the coo l, dense, post-shock layers not availabl e in other 
SNe (e.g.. lPozzo etaDl200lt ISmith et al.l[2009cT ). and (2) 
the dense CSM produced by pre-SN mass loss can con- 
tain dust that is illuminated by the superno va radiation 
(e.g.. lSlnitbetld1[2^cllMiller et al.ll2010aD . 

The origin and heating mechanism of the dust, how- 
ever, can be ambiguous. Disentangling the various dust 
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models offers important clues regarding the supernova's 
circumstellar environment, explosion mech anism, and 
even progenitor system (e.g.. lFox et al.lr2"01Q[ ). For exam- 
ple, if shock heat ed, the dust temperature c orrelates with 
the ga s den sity (iDraine fc Salpeterl [19791: IDraind IT98TI: 
iDwekl 119871 : iDwek et al.ll2008f) . which can be used to 
trace the progenitor's m ass-loss history (e.g. . ISmith et all 
I2009d : IFox et al.ll2010h . Alternatively, the discovery of 
significant amounts of newly formed SN dust would pro- 
vide the much sou ght evidence necessa r y to c onfirm SN 
dust models (e.g.. iNozawa et al J [20031 120081 ) . At this 
point, the observed dust yields all tend to be 2-3 or- 
ders of magnitude smalle r than predicted by the models 
(e.g.. iKozasa et al.|[200l iMeikle etail I2011L and refer- 
ences therein), although far-IR observations of SN 1987A 
suggest larger reservoi rs of dust may be hidd en at colder 
temperatures (~20 K. IMatsuura et aT1l2011l ). 

Mid-IR observations, which span the peak of the ther- 
mal emission from warm dust, offer the most useful con- 
straints on the dust properties. Yet obs ervations a t these 
wavelengths remain sparse. In fact, the IFox et all ([20101 ) 
Spitzer Space Telescope IRS spectrum of SN 2005ip is 
the only mid-IR spectrum of a SN Iln to date. Es- 
timated to h ave a rate of no more th an 10 yr _1 out 
to 150 Mpc (jDahlen fc Franssonl [1999T ) . the number of 
nearby SNe Iln is small, and SNe Iln at greater distances 
are not easy targets. 

The few observed dust-emitting SNe Iln, however, are 
anomalously luminous and linger for many years at mid- 
IR wavelengths, thereby making remnant archeology pos- 
sible. While 79 SNe Iln within a distance of 250 Mpc 
were discovered during the years 1999-2008, only 11 
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Table 1 

nth Observed Late-Time IR Excess 



SN 



Subclass 



Spitzer Observations?/ 
Detections? 



Reference 



1979c 

1980k 

1982e 

19821 

1982r 

19851 

1986J 

1987A 

1993J 

1994Y 

1995G 

1995N 

1997ab 

1998S 

1999Z 

1999bw 

1999el 

1999em 

2002bu 

2002hh 

2002ic 

2003gd 

2004dj 

2004et 

2005af 

2005kd 

2005ip 

2006gy 

2006jc 

2006tf 

2007it 

2007od 

2007rt 

2008iy 



11 
II 

? 

II 
I 

11 
lln 
II 
lib 

? 

lln 
lln 
lln 
lln 
lln 
lln 
lln 
IIP 
lln 
IIP 
Ia/IIn 
IIP 
IIP 
IIP 
IIP 
lln 

IIll 

lln 
lb 
lln 
IIP 
IIP 
lln 
lln 



No/No 
No/No 
No/No 
No/No 
No/No 
No/No 
No/No 
Yes/ Yes 
No/No 
No/No 
No/No 
No/No 
No/No 
No/No 
No/No 
Yes/ Yes 
No/No 
No/No 
Yes/ Yes 
Yes/ Yes 
No/No 
Yes/Yes 
Yes/ Yes 
Yes /Yes 
No/No 
No/No 
No/No 
No/No 
No/No 
No/No 
Yes/ Yes 
No/No 
No/No 
No/No 



Merrill (19801 
Telesco et aT"Tl9811 
Graham et al. (1983) ; Graham"*" Meikle (1986' 
Graham fc Meikle (1986^ 
Graham fc Meikle (Tf*3*5l 
Elias et al. (1986) ; Fesen7l998') 
Weiler^Tal. ( 1990O*eibundgut et al. (19911 
Many Authors 
Lewis et al. (1994) ; Matheson et al. (20001 
Garnavich""3t""d" 7**j9** " — 
Pastorello et al. (2002*) 
Gerardv et al. (2002) 
Gerardv et al. (2002) 
Gerardv et al. ( 2002 ) "Pozzo et al. (20041 
Gerardv et al. ( 200*2* "" 
Sugerman et al. (5*J04l 
Gerardv et al. (2002) *Carlo et al. (20021 
Elmhamdi et al. (2003) 
Thompson et al. (20l59l 
Meikle et al. (2o75cT) 
Kotak et al. (2CJ0T) 
Sugerman et al. f2006**~Meikle et al. (20071 
Kotak et al. (2005 , 2006); Szalai et al. (20111 
Kotak et alH**%9*"*Maguire et al. (20T**T 
Kotak et al. (2006* "" 
Tsvetkov (20087 
Fox et al. (2009) ; Smith et al. (2009c 1 
Smith et al. (2009aT ' 
Smith et a~ (2008b) 
Smith et al. ^008a) 
Andrews et al. (20111 
Andrews et al. (2010) 
Trundle et al (2009) 
Miller" 3t al. (2010b) 



of these positions have archived Spitzer IRAC data 
collected (either intentionally or serendipitously) more 
than a few months post-discovery. Five of these eleven 
events (^45%) show evidence for late-time IR emission, 
in some cases up to six years post-discovery: SN e 1999bw 
( Sugerman eFatl l200i IThompson eTaLI [2001) 1999 el 
( Carlo et all [2002ll. 2002 bu ([Thompson et all [2009lf| . 
20031 o (jMeikle et al.ll2005l) . and 2005ip (|Fox et al.ll2009l 
[20101) . 

We therefore executed a warm Spitzer follow-up sur- 
vey of the remaining, 68 Type lln events discovered in 
the years 1999-2008Q This paper presents the results of 
the survey, including the positive detection of late-time 
emission from ten targets, which more than doubles the 
database of existing mid-IR observations of SNe lln. For 
some targets, optical and near-IR photometry and spec- 
troscopy also exist. In JjHwe summarize the observations 
and data-reduction techniques. Spitzer photometry con- 
strains the dust mass, temperature, and, thereby, lumi- 
nosity. Section EH implements the methods presented by 
IFox et all ^OloT to explore the origin and heating mech- 
anism of these components; in most cases, results are 
compared to those of SN 2005ip. Section 2] discusses the 
overall trends, circumstcllar model, and progenitor sys- 
tem, while fj5] summarizes the results and discusses future 



work. 



2. OBSERVATIONS 



6 Note that ISmith et all ||20TO|1 suggest that SNe 1999bw and 
2002bu may be SN impostors instead of genuine SNc. 

7 We also included the well-studied Type lln SN 1997ab, for a 
total of 69 SNe. 



2.1. Mid-Infrared Warm Spitzer IRAC Survey and 
Photometry 

Warm Spitzer /IRAC (jFazio et al.ll2004[) surveyed the 
positions of the 69 SNe lln in the sample. Although 
the Spitzer warm mission is limited to IRAC bands 1 
(3.6 fjLm) and 2 (4.5 /xm), ground-based observations can- 
not compete because they have limited sensitivity past 
2 ivco. given the rapidly rising thermal background. Fur- 
thermore, these two bands span the peak of the black- 
body emission from warm dust (500 < < 1000 K), 
providing potential constraints on the dust temperature 
and mass (and thus the luminosity). Assumin g peak lu- 
mino sities comparable to that of SN 2005ip (|Fox et al.l 
I2010D for each SN, we developed Astronomical Observa- 
tion Requests (AORs) that would yield signal-to-noise 
ratio (SNR) > 25 at 4.5 /im to ensure the most-reliable 
extraction of the temperature and mass from the fluxes 
in the two bands. All integration times were rounded up 
to a multiple of 5 min to produce a simple and straight- 
forward AOR. Table lists the observational details. 

The LeoparcQ software package provided access to the 
Post Basic Calibrated Data (pbed), which are already 
fully coadded and calibrated. Supernova detections were 
initially identified with photometry routines (SNR > 3) 



8 Leopard can be downloaded 

http://ssc.spitzcr.caltcch.edu/propkit/spot/ . 



from 
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Table 2 

Spitzer Survey Targets 



SN 


JD 


Epoch 


a (h:m:s) 


<5 (°:':") 


Distance 


'int 




(-2,450,000) 


(days) 


(J2000) 


(J2000) 


(Mpc) 


(s) 


1997ab 


5140 


5018 


09:51:00.40 


+20:04:24.0 


50 


300 


1999Z 


5065 


4104 


10:22:37.23 


+27:21:19.8 


210 


1800 


1999eb 


5170 


3612 


01:43:45.45 


+04:13:25.9 


75 


300 


2000P 


5322 


3450 


13:07:10.53 


-28:14:02.5 


31 


300 


2000cl 


5063 


3631 


10:37:16.07 


-41:37:47.8 


38 


300 


2000ct 


5066 


3335 


17:01:03.64 


+33:28:45.0 


124 


300 


2000eo 


5068 


3247 


03:09:08.17 


-10:17:55.3 


44 


300 


2000ev 


5174 


3264 


06:47:52.00 


+84:10:02.2 


61 


300 


20011 


5322 


3180 


03:43:57.28 


+39:17:39.4 


68 


300 


2001ey 


5074 


3002 


22:26:30.77 


-06:23:28.6 


105 


300 


2001fa 


5065 


2872 


01:48:22.22 


+11:31:34.4 


71 


300 


2001ir 


5209 


2937 


08:36:28.12 


-11:50:03.5 


81 


300 


2002A 


5075 


2864 


07:22:36.14 


+71:35:41.5 


40 


300 


2002bj 


5209 


2784 


05:11:46.41 


-15:08:10.8 


50 


300 


2002bv 


5146 


2796 


07:50:01.13 


+30:01:32.2 


115 


300 


2002cb 


5186 


2950 


13:04:23.97 


+47:35:53.2 


123 


300 


2002ca 


5118 


2584 


02:08:25.08 


+ 14:20:52.8 


61 


300 


2002fj 


5108 


2792 


08:40:45.10 


-04:07:38.5 


60 


300 


2003as 


5078 


2455 


05:28:45.81 


+49:52:59.1 


95 


300 


2003ci 


5095 


2296 


16:28:40.30 


+ 12:46:05.8 


110 


300 


2003hy 


5084 


2290 


21:54:22.72 


+ 15:09:37.8 


101 


300 


2003hy 


5084 


2290 


21:54:22.72 


+ 15:09:37.8 


101 


300 


2003jh 


5202 


2166 


04:16:37.30 


-12:23:32.9 


122 


300 


2004F 


5151 


2067 


03:17:53.80 


-07:17:43.0 


73 


300 


2004ay 


5065 


1977 


18:28:57.57 


+51:38:55.7 


133 


300 


2004ec 


5322 


1825 


17:08:04.96 


+26:22:44.4 


171 


900 


2004gd 


5113 


1835 


07:09:11.71 


+20:36:10.6 


72 


300 


2005R 


5075 


1938 


11:14:56.81 


+33:49:35.7 


151 


600 


2005aq 


5135 


1671 


04:31:38.82 


-04:35:06.8 


56 


300 


2005av 


5066 


1657 


20:44:37.58 


-68:45:10.6 


13 


300 


2005bx 


5065 


1711 


13:50:24.95 


+68:33:19.4 


132 


300 


2005cl 


5084 


1635 


21:02:02.35 


-06:17:35.7 


107 


300 


2005cp 


5066 


1523 


23:59:30.88 


+18:12:09.6 


91 


300 


2005cy 


5066 


1498 


18:26:49.11 


+51:08:30.4 


138 


300 


2005db 


5111 


1495 


00:41:26.79 


+25:29:51.6 


62 


300 


2005dr 


5108 


1476 


05:02:17.18 


+07:38:21.4 


157 


600 


2005ds 


5322 


1470 


17:57:39.89 


+27:50:18.4 


97 


300 


2005cp 


5151 


1421 


16:56:35.49 


+58:01:30.2 


121 


300 


2005gl 


5072 


1417 


00:49:50.02 


+32:16:56.8 


65 


300 


2005gn 


5065 


1479 


05:48:49.07 


-24:22:45.5 


165 


600 


20051y 


5108 


1349 


01:23:28.84 


+30:46:45.7 


143 


600 


2005ma 


5187 


1384 


04:49:53.91 


-10:45:23.4 


62 


300 


2006aa 


5075 


1547 


11:53:19.89 


+20:45:18.2 


85 


300 


2006am 


5139 


1276 


14:27:37.24 


+41:15:35.4 


37 


300 


2006bo 


5322 


1320 


20:30:41.90 


+09:11:40.8 


63 


300 


2006bv 


5077 


1348 


12:41:01.55 


+63:31:11.6 


35 


300 


2006cu 


5322 


1181 


14:47:43.31 


+09:39:33.9 


119 


300 


2006cy 


5159 


1426 


13:08:01.23 


+26:06:59.0 


150 


600 


2006dc 


5118 


1179 


16:16:04.04 


-22:37:16.7 


107 


300 


2006ch 


5140 


1127 


02:50:11.95 


+47:03:20.5 


111 


300 


2006fp 


5200 


1082 


22:45:41.13 


+73:09:47.8 


21 


300 


2006gy 


5073 


1111 


03:17:27.06 


+41:24:19.5 


76 


300 


2006jd 


5201 


1149 


08:02:07.43 


+00:48:31.5 


77 


300 


2006qq 


5108 


1048 


05:19:50.43 


-20:58:06.4 


119 


300 


2007K 


5081 


1069 


09:14:26.21 


+36:06:52.4 


90 


300 


2007ak 


5322 


965 


05:20:40.75 


+08:48:16.0 


66 


300 


200 1 bb 


5062 


953 


0/:01:0/.46 


+ol:lo:o/.o 


86 


300 


2007cm 


5066 


964 


12:42:45.18 


+55:08:57.1 


66 


300 


2007pk 


5322 


660 


01:31:47.07 


+33:36:54.1 


69 


300 


2007rt 


5203 


780 


11:02:34.29 


+50:34:58.5 


93 


300 


2008B 


5066 


597 


15:02:43.65 


+23:20:07.8 


79 


300 


2008J 


5135 


593 


02:34:24.20 


-10:50:38.5 


66 


300 


2008aj 


5322 


812 


13:33:06.33 


+33:08:59.0 


102 


300 


2008bc 


5066 


636 


14:27:48.71 


+69:41:50.8 


123 


300 


2008cg 


5135 


474 


15:54:15.15 


+ 10:58:25.0 


151 


600 


2008cn 


5065 


386 


00:55:13.56 


+35:26:26.2 


151 


600 


2008fm 


5169 


351 


23:49:03.51 


+26:47:39.3 


157 


600 


2008gm 


5199 


301 


23:14:12.39 


-02:46:52.4 


49 


300 


2008ip 


5088 


490 


12:57:50.20 


+36:22:33.5 


63 


300 



All distances are derived from the host-galaxy redshift assuming Ho 
= 72 km s" 1 Mpc -1 . 
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Figure 1. False-color 3.6 (im and 4.5 (im combined Spitzer /IRAC 
images of the 10 SNe Iln discovered by our survey (P60122) to 
exhibit late-time mid-IR emission. 

and later confirmed with a visual inspection. Of the 69 
SN positions surveyed, f SNe exhibit late-time emission 
(see Figure [T]). 

The background flux in most of the SN host galaxies 
is bright and exhibits rapid spatial variations, making 
faint detections difficult for the DADPHDT point-spread 
function (PSF) photometry package in IRAF0 and lim- 

9 IRAF: the Image Reduction and Analysis Facility is distributed 
by the National Optical Astronomy Observatory, which is oper- 



SN 2005cp 
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Figure 2. Illustration of the aperture-photometry routine im- 
plemented instead of DAQPHOT and APPHDT (which utilizes both an 
inner aperture and outer annulus). A set of single apertures, iden- 
tified by the green circles, extracts the flux from both the SN and 
a user-selected background. This technique allows the user to vi- 
sually identify and extract local background representative of the 
underlying SN flux. Nearby stellar and nuclear flux can be ignored. 

iting the accuracy of sky-brightness measurements using 
the annuli in the IRAF APPHOT (aperture photometry) 
package. Although template subtraction is a commonly 
used technique to minimize photometric confusion from 
the underlying galaxy, no pre-SN observations exist. 

Instead, photometry was performed by using a set of 
single apertures, with a radius defined by a fixed mul- 
tiple of the PSF full width at half-maximum intensity 
(FWHM), to estimate both the SN and average back- 
ground flux. This technique allowed the user to visually 
identify only local background associated with the SN, 
as opposed to the annuli imposed by APPHOT. Figure [5] il- 
lustrates this technique in the case of SN 2005cp, where 
the circles identify the apertures used to extract both 
the SN and background fluxes. Nearby stellar and nu- 
clear flux, which would normally be included in the sky 
annuli, do not significantly contribute to the SN flux and 
are ignored by this technique. 

Pixel fluxes were converted from mJy sr _1 to mJy ac- 
cording to the IRAC Data Handbook version 3, which 
discusses the pixel size and aperture correction in de- 
tail. The DADPHDT and APPHOT results for the brighter 
SNe agreed within the uncertainties. Because DA0PH0T 
did not succeed in detecting the fainter events, however, 
we only include results from APPHOT. These fluxes are 
listed in Table [3] and plotted in Figure [3] Upper limits 
for nondetections (~0.015 mJy at 4.5 /im) were set by 
the point-source sensitivity in Table 2.10 of the IRAC 
Instrument Handbook version 2. In some cases, coinci- 
dent near-IR data from several other instruments exists 
(see below). 

2.2. Spitzer Survey Supernova Detections 

SN 2005cp was discovered in UGC 12886 on 2005 
June 20.48 (UT d ates are u s ed th roughout this paper; 
R w 17. 4 mag) by | Lee et al.l (|2005D and spectrally clas- 
sified by ^^m^AT^M) ■ The line profiles consist 

ated by the Association of Universities for Research in Astronomy 
(AURA) under cooperative agreement with the National Science 
Foundation (NSF). 
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Table 3 

Spitzer Survey Detections 1 



SN 


Epoch 


* 

3.6 fim 


4.5 fim 




(days) 


(mjy) 


(mjy) 


2005cp 


1523 


0.25(0.07) 


0.31(0.07) 


2005gn 


1479 


0.11(0.04) 


0.12(0.04) 


2005ip 


948 


5.76(0.14) 


7.54(0.22) 


2006jd 


1149 


2.00(0.17) 


2.71(0.20) 


2006qq 


1048 


0.14(0.07) 


0.24(0.07) 


2007rt 


780 


1.60(0.15) 


1.96(0.17) 


2008J 


593 


2.52(0.18) 


2.71(0.19) 


2008cg 


474 


0.30(0.07) 


0.35(0.08) 


2008en 


386 


0.30(0.07) 


0.34(0.07) 


2008gm 


301 


0.05(0.04) 


0.03(0.03) 


2008ip 


438 


0.05(0.03) 


0.07(0.03) 



lcr uncertainties are given in parentheses. 
1 Upper limits for nondetections (^0.015 
mjy at 4.5 fj,m) were set by the point-source 
sensitivity in Table 2.10 of the IRAC Instru- 
ment Handbook version 2. 
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Figure 3. Photometry of the SNe detected by Spitzer /IRAC 
bands 1 (3.6 fim) and 2 (4.5 p.m), as well as near-lR photometry ob- 
tained by Magellan/PANIC. Overplotted are the APO/TripleSpec 
near-IR spectra of SN 2006jd and SN 2005ip. 



of both a narrow (FWHM w 400 km s _1 ) and interme- 
diate (FWHM w 4000 km s _1 ) component. iKiewe et all 
()2010[) provide additional early-time optical photometry 
and spectra, as well as an associated mass-loss rate. 

SN 2005gn was discovered in E SQ 488-G30 on 2005 
Oct ober 13.53 (R j 18 mag) by ILuckas et all (|2005h 
and IPrasad fc Lil (|2005 ), and s p ectra lly classified by 
iBlanc et all (|2005f ) and lMarsdenl (|2005l ). The line pro- 



files consist of narrow, intermediate, and broad features, 
although no line widths were reported. 

SN 2006jd was discovered in U GC 4179 on 2006 Oc - 
tob er 12.54 (R « 17.2 mag) by lltagaki et all ([20061 ) 
and IPrasad feU (f2006h . It was classified as a SN Iln, 
but no line widths were reported. Ilmmler et al.l ()2007l ) 
found an associated X-ray luminosity of Lx — 2.5 x 
10 41 erg s _1 and a visual magnitude V = 19.1 on 2007 
November 16.73. Swi/Z/XRT measured an X-ray lumi- 
nosity Lx = 2.5 x 10 
|2007|) . and lChandra 



41 



erg s 1 on day 4 00 (Ilmmler et al.l 



SN 2006qq was discovered in ESO 553-G36 on 2006 
October 27.40 (R » 17.0, M v = -18.4 mag) by 
IPrasad et al.l (|2006[ ) and classified by iSilverman et al.l 
(|2006| ). The line profiles consist of narrow (FWHM w 
900 km s _1 ) emission upon a relatively featureless con- 
tinuum. 

SN 2007rt was discovered in UGC 6109 on 2007 
November 24.54 (R « 16.8 mag) by ()2007l ). It was 
spectrally classi fied as a SN Iln about 2- 3 months past 
maximum light (jBlondin fe Berlindll2007i) . The line pro- 
files consist of both an intermediate (FWHM w 2500 
km s" 1 ) and broad (FWHM w 10,000 km s" 1 ) compo- 
nent. A moderate-resolution spectrum reveals an addi- 
tional narrow Ha P-Cygni profile (FWHM » 100 km s" 1 ) 
and evidence for new dust formation dTrundle et al.l 
[2001 . On dav 48. IChandra fc Soderberg! (|2008br i did not 
detect any radio emission and assigned an upper limit of 
< 9 ptJy 

SN 2008J was discovered in MCG-02-7-33 on 2008 
February 15.19 (R . » 15.4 mag) bvlThrasher et al.l (|2008l ) 
and classified by dStritzinger et al.l l2008f ) . The line pro- 
files consist of both a narrow and intermediate compo- 
nent, although the line wi^h^were not listed. 

SN 2008cg was discoveredin FGC 1965 on 2008 May 
5.39 (R k, 16.7) bvlDrake et al.l (|2008f ) and classified by 
(|Blondin fc Hernan dez 2008). The line profiles consist of 
both a narrow (FWHM w 300 km s" 1 ) and intermediate 
(FWHM w 1500 km s" 1 ) compo nent. ffw/VUVOT ' mea- 
sured V = 17.9 mag on dav 58 (Ilmmler et al.ll2008f ). On 
day 57, however. I Chandra fc Soderberg! ( 2008a! ) did not 



flux density of 238 fiJy on day 405. 



Soderberg ( 2007) reported a radio 



detect any radio emission and assigned an upper limit of 

< 7 /uJy. 

SN 2008en was discover ed in UGC 5 64 on 2008 Au- 
gust 3.06 (R « 18.2 mag) bvlBoled (|2008| ) and spectrally 
classified by iSteele et al.l (|2008l ). These spectra reveal 
broad (FWHM « 3700 km s _1 ) Balmer and He I emis- 
sion lines. 

SN 2008gm was discovered in NGC 7530 on 2008 
October 22.09 ( R » 17.0 mag ) bv iPignata et al.1 (poM ) 
and classified by iPrieto et al.l (|2008l ). The line profiles 
consist of both a narrow (FWHM w 300 km s _1 ) and 
intermediate (FWHM w 1600 km s -1 ) compon ent. Two 
days following the discovery, iSoderbergl ()2008l ) did not 
detect any radio emission and assigned an upper limit of 

< 0.08 mjy (3a). 

SN 2008ip was discovered in N GC 4846 on 2008 De - 
cember 31.75 (R g 15.7 mag) bv [Nakano et al.l (|2009T > 
and classified by iChallis fc Calkins! (I2009D. No spec tral 
line widths were given. iChandra fc Soderberg] (|2009f ) re- 
port no radio counterpart on day 5 and assign a 3<7 upper 
limit of 0.096 mjy (8.46 GHz). 

2.3. Near-Infrared Photometry and Spectroscopy 

Nearly coincident near-IR photometric data were col- 
lected for several targets with th e Persson Auxiliary 
Nasmyth Infrared Camera (PANIC; Partini et al.ll2004l) 
mounted on the 6.5-m Magellan (Walter Baade) Tele- 
scope at the Las Campanas Observatory, Chile. PANIC 
obtained J- and K s -band photometry of SNe 2008gm 
(363 days post-discovery), 2008cg (659 days), 2006jd 
(1230 days), and 2005cp (1582 days). Standard near- 
IR data-reduction techniques were implemented. The 
relatively small galaxy sizes allow data frames to dou- 
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Table 4 

Summary of Optical Spectra 



Table 5 

Summary of Keck/LRIS Optical Photometry 



SN 


Epoch 

(days) 


Instrument 


Blue 
Resolution 


Red 
Resolution 


Int 

00 


2005cp 


10 


Lick/Kast 


7.7 


10.4 


900 


2006jd 


395 
564 


Keck/LRIS 
Keck/LRIS 


6.5 
6.5 


7 

5.8 


900 
600 


2006qq 


1 

20 


Lick/Kast 
Lick/Kast 


6 

6.7 


10.9 
10.0 


1800 
2100 



SN 





32 


Lick/Kast 


4.9 


11.9 


900 




199 


Lick/Kast 


6 


11 


900 




224 


Lick/Kast 


6.7 


12.3 


900 


2008J 


236 


Lick/Kast 


7.1 


12.3 


1200 




251 


Lick/Kast 


4.6 


11.0 


1200 




266 


Lick/Kast 


1,1 


11.0 


1200 




313 


Lick/Kast 


4.8 


5 


1200 




3 


Lick/Kast 


5.3 


11.6 


1800 




10 


Lick/Kast 


5.5 


11.7 


1800 


2008cg 


55 


Lick/Kast 


5.9 


10.8 


1500 




63 


Lick/Kast 


6 


11.3 


1500 




114 


Keck/LRIS 


4.5 


7 


454 


2008en 


5 


Lick/Kast 


6.7 


12.5 


1800 


2008gm 


5 


Keck/LRIS 


6.5 


7 


400 



ble as sky exposures, provided a sufficiently large dither 
avoids an overlapping galaxy from one exposure to the 
next. Photometry was performed with the DAOPHOT 
PSF-fitting techniques. The results are plotted in Fig- 
ure [3] adjacent to the Spitzer photometry and, together, 
are treated as a single epoch given the near simultane- 
ity of the observations relative to the overall epoch of 
observation. 

TripleSpcc, a 0.9-2.5 /im, resolution R = 3000 spec- 
trograph operating at Apache Point O bservatory in NM 
(| Wilson et all 120041: IHerter et all 120081 ). was used to ob- 
tain a spectrum of SN 2006jd on day 1250 post-discovery 
and a spectrum of SN 2 005ip on day 86 2 post-discovery 
(originally published bv lFox et al.ll2010l ). Forty minutes 
of on-source integration consisted of 8 independent 5 min 
exposures nodding between 2 different slit positions. We 
extracted the spe ctra with a modified version of the IDL- 
based SpexTool (|Cushing et al.l |2"004[) . Any underlying 
galactic and sky emission are approximated in SpexTool 
by a polynomial fit and subtracted from the SN. The re- 
sults are again plotted in Figure [3] and, combined with 
the Spitzer data, treated as a single epoch. 

2.4. Optical Photometry and Spectroscopy 

Table [?] summarizes the optical spectra obtained with 
various instruments. Data were obtained with the 
dual-arm Low R esolution Imaging Spectrometer (LRIS; 
lOke et all I1995D mounted on the 10-m Keck I tele- 
scope, as well as wi th the Kast double spectrograph 
(|Miller fc Stond 119931 ) mounted on the Shane 3-m tele- 
scope at Lick Observatory. Keck/LRIS spectra were ob- 
tained using the 600/4000 or 400/3400 grisms on the blue 
side and the 400/8500 grating on the red side, along with 
a 1" wide slit. This resulted in a wavelength coverage 
of 3200-9200 A and a typical resolution of 5-7 A. Al- 
most every Lick/Kast spectrum was obtained using the 
600/4310 grism on the blue side and the 300/7500 grating 
on the red side, along with a 2" wide slit. This resulted in 
a wavelength coverage of 3300-10,400 A and resolutions 
of about 5-7 A and 10-12 A on the blue and red sides, 
respectively Most observations had the slit aligned along 



Epoch 
(days) 



B 
(mag) 



R 

(mag) 



(mag) 



Luminosity 
(log (L opt /L )) 



2006jd 1609 
2007rt 1223 



25.0±0.5 

>25 



20.3±0.5 
>24 



22.5±0.5 
>24 



7.6±0.4 
<6.1 




Figure 4. False-color BRI combined Keck/LRIS images of SN 
2006jd on day 1609 post-detection. The SN remains bright at op- 
tical wavelengths nearly 4.5 yr after the explosion. 



the parallactic an gle to minimize differential light losses 
(lFiliDDenkolll982D . 

Th e spectra were red uced using standard techniques 
(e.g.. iFolev eFaII[200l . Routine CCD processing and 
spectrum extraction were completed with IRAF, a nd the 
data w ere extracted with the optimal algorithm of lHornel 
(|1986| ) . We obtained the wavelength scale from low-order 
polynomial fits to calibration-lamp spectra. Small wave- 
length shifts were then applied to the data after cross- 
correlating a template sky to an extracted night-sky spec- 
trum. Using our own IDL routines, we fit a spectropho- 
tomctric standard-star spectrum to the data in order to 
flux c alibrate the SN and to remove telluric abs orption 
lines (|Wade fc Hornd[l98l iMatheson et al.ll200Q[ ). 

Kcck/LRIS also obtained BVRI photometry of SNc 
2006jd (1609 days) and 2007rt (1223 days) on 2011 
March 9. Observations of each target consisted of two 
dithered 120-s exposures in each band. Again, the 
data were reduced using standard techniques and routine 
CCD processing with IRAF. Table [5] lists the photomet- 
ric results. Only SN 2006jd was detected (see Figured]), 
but upper limits for SN 2007rt are also provided. 

2.5. Dust Composition, Temperature, and Mass 

Assuming only thermal emission, the combined near- 
IR and mid-IR spectra provide a strong constraint on the 
dust temperature and mass, and thus on the IR luminos- 
ity. Assuming optically thin dust with mass M d , with a 
particle radius a, at a distance d from the observer, ther- 
mally emitting at a single equi librium temperat ure 
the flux can be written as (e.g.. lHildebrandlll983f ) 



M d B v (T d ) Kv (a) 
d 2 
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X (jim) X (fim) 

(a) Graphite (b) Silicate 

Figure 5. Best fits of Equation [T] (both graphite and silicate compositions) to the combined near-IR and mid-IR data. A warm (~ 400- 
800 K) dust component fits most observations well. SNe 2006jd and 2005ip (those with near-IR spectra) are better fit by a multi-component 
model, where a second, hotter component dominates the near-IR flux. The most obvious feature differentiating between the graphite and 
silicate models is the 9 fim silicate emission feature. The warm IRAC data do not extend beyond ~ 5 (im, limiting any distinction between 
the two compositions. 



Table 6 

Spitzer Fitting Parameters 







Graphites 




Silicates 






Blackbody 




a (fim) 


T d (K) 


M d (M ) 


L d (L ) 


T d (K) M d (Mq) 


L d (Lq) 


T d (K) 


M d (M ) 


L d (L ) 


SN 2005cp 


O.Of 


575 


7.39e-3 


5.65e+7 


648 1.51e-2 


2.06c+8 


710 


9.09e-6 


6.86e+7 


0.10 


554 


6.81e-3 


5.63e+7 


643 1.55e-2 


2.07e+8 


710 


9.08e-5 


6.85e+7 


0.30 


580 


1.79e-3 


5.15e+7 


614 1.78e-2 


2.13e+8 


710 


2.70e-4 


6.79c+7 


SN 2005gn 


O.Of 


776 


2.06e-3 


8.39e+7 


991 3.12e-3 


1.72e+8 


1220 


1.52e-6 


1.06e+8 


O.fO 


737 


1.92e-3 


8.36e+7 


981 3.18e-3 


1.73e+8 


1220 


1.51e-5 


1.06e+8 


0.30 


735 


6.75e-4 


7.52e+7 


917 3.67e-3 


1.75e+8 


1220 


4.50e-5 


1.05e+8 


SJN 2005ip 


O.Of 


483 


4.73e-2 


1.36e+8 


525 9.12e-2 


6.23e+8 


567 


5.f2e-5 


1.49e+8 


O.fO 


468 


4.34e-2 


1.35e+8 


522 9.30e-2 


6.26e+8 


567 


5.ffe-4 


1.48e+8 


0.30 


469 


1.46e-2 


1.19e+8 


503 1.07c-l 


6.50c+8 


567 


f.52e-3 


1.47e+8 


SN 2006jd 


O.Of 


602 


3.21e-2 


3.18e+8 


733 4.90e-2 


9.94e+8 


836 


2.70e-5 


4.03e+8 


O.fO 


579 


2.96e-2 


3.17e+8 


728 5.00e-2 


1.00e+9 


836 


2.70e-4 


4.02e+8 


0.30 


561 


1.13e-2 


2.65e+8 


692 5.76e-2 


1.02e+9 


836 


8.0fe-4 


3.98e+8 


SN 2006qq 


O.Of 


518 


1.71e-2 


7.30e+7 


601 2.75e-2 


2.94e+8 


672 


f.45e-5 


8.68e+7 


O.fO 


501 


1.58e-2 


7.27e+7 


598 2.79e-2 


2.95e+8 


672 


f.45e-4 


8.66e+7 


0.30 


500 


5.54e-3 


6.66e+7 


575 3.16e-2 


3.03c+8 


672 


4.30e-4 


8.58e+7 


SN 2007rt 


O.Of 


653 


2.39c-2 


3.71e+8 


793 3.76e-2 


9.84e+8 


926 


1.92e-5 


4.37e+8 


O.fO 


625 


2.22e-2 


3.70e+8 


787 3.83e-2 


9.87e+8 


926 


1.92e-4 


4.36e+8 


0.30 


624 


7.79e-3 


3.41e+8 


747 4.37e-2 


1.00e+9 


926 


5.70e-4 


4.32e+8 


SN 2008J 


O.Of 


739 


9.41e-3 


2.92e+8 


930 1.44e-2 


6.37e+8 


1125 


7.13e-6 


3.58e+8 


O.fO 


704 


8.73c-3 


2.91e+8 


921 1.47e-2 


6.40e+8 


1125 


7.12e-5 


3.58e+8 


0.30 


702 


3.07e-3 


2.64e+8 


865 1.69e-2 


6.48e+8 


1125 


2.12e-4 


3.54e+8 


SN 2008cg 


O.Of 


629 


1.43e-2 


1.80e+8 


716 2.93e-2 


5.52e+8 


975 


7.70e-6 


2.17e+8 


O.fO 


604 


1.32e-2 


1.79e+8 


710 3.02e-2 


5.57e+8 


975 


7.69e-5 


2.16e+8 


0.30 


635 


3.43e-3 


1.67e+8 


675 3.48e-2 


5.70e+8 


975 


2.28e-4 


2.14e+8 


SN 2008en 


O.Of 


714 


7.23e-3 


1.84e+8 


888 l.lle-2 


4.23e+8 


1062 


5.58e-6 


2.23e+8 


O.fO 


681 


6.70e-3 


1.84e+8 


881 1.14e-2 


4.25e+8 


1062 


5.57e-5 


2.22e+8 


0.30 


679 


2.36e-3 


1.68e+8 


829 1.30e-2 


4.31c+8 


1062 


1.66c-4 


2.20e+8 


SJN 2008gm 


O.Of 


757 


7.00e-5 


2.47e+6 


890 1.40e-4 


5.36e+6 


1013 


8.31e-8 


2.73e+6 


O.fO 


721 


6.47e-5 


2.46e+6 


880 1.45e-4 


5.39e+6 


1013 


8.29e-7 


2.73e+6 


0.30 


766 


1.69e-5 


2.36e+6 


827 1.68e-4 


5.48e+6 


1013 


2.46e-6 


2.70e+6 


SJN 2008ip 


O.Of 


634 


4.52e-4 


5.93e+6 


765 7.12e-4 


1.66e+7 


886 


3.66e-7 


6.96e+6 


O.fO 


608 


4.18e-4 


5.91e+6 


759 7.25e-4 


1.66e+7 


886 


3.66e-6 


6.95e+6 


0.30 


606 


1.47e-4 


5.46e+6 


721 8.27e-4 


1.69e+7 


886 


1.09e-5 


6.88e+6 



8 



where _B„(Td) is the Planck blackbody function and K v {a) 
is the dust mass absorption coefficient. 

For simple dust populations of a single size composed 
entir ely of either silicat e or graphite, the IDL MPFIT func- 
tion (|Markwardtll2009fl finds the best fit (see Figure 
of Equation Q] by varying and T<j to minimize the 
value of y 2 . The abs o rption coefficients, n, are given in 
Figure 4 of iFox et all (|2010D . The combined photometry 
and near-IR spectra of SNe 2006jd and 2005ip are best 
fit by a multi-component model. As opposed to a single 
component with a continuous temperature distribution, a 
second, hotter component at a single temperature dom- 
inates the near-IR flux, while also contributing to the 
3.6 /jm flux. For the remaining SNe, the small num- 
ber of photometric data points limits our fits to a single 
component. Since a second, hotter component may be 
contributing to the shorter wavelength flux, we treat the 
derived dust temperatures for these SNe as upper limits 
(see 33.31 below for significance). 

Table [6] lists the best-fit parameters (i.e., temperature 
and mass, as well as IR luminosity) for several grain 
sizes. Overall, grain size has little consequence on the 
goodness of fit or the resulting parameters. The only 
size- dependent param eter in Equation[T]is k, but Figure 
4 of IFox et al.l (|2010f ) shows that the dust-opacity coef- 
ficient for both graphite and silicate is independent of 
grain radius at IR wavelengths (>1 /an) for grain sizes 
<1 ^tm, which is typical for most grains. A grain size of 
a = 0.1 fim is therefore assumed throughout the rest of 
this paper, including the results plotted in Figure O 

In general, the derived dust temperatures, masses, and 
luminosities are all higher if we assume silicate grain com- 
position. The limited number of data points (2 or 3 in 
most cases), however, limits the best-fitting function in 
Figure [5] from distinguishing between the two composi- 
tions. The most obvious differentiating feature is the 9 
fim silicate emission band, but the warm IRAC data do 
not extend beyond ~5 /im. Both compositions are there- 
fore considered throughout the rest of this paper. 

Figure |6] plots the temperature, mass, and IR luminos- 
ity for a = 0.1 /im grains for both the graphite and sili- 
cate fits, as well as the near- IR evolution of the Type Iln 
SNe 2005ip (IFox et al.ll20"09l) and 1995N (IGerardv et al.l 
120021) . The luminosity curve plateaus to a point between 
1000 and 2000 days (^3-5 yr), at which time it seems to 
slowly fade with an e-fold time of several hundred days. 
Although many of the points are upper limits, the trend 
appears to closely follow the observed near-IR decline of 
SN 1995N. Furthermore, at least half of the SNe have 
derived dust masses > 10 -2 M© (and most of the 
rest have masses Afd > 10~ 3 M Q ). If newly formed, 
these masses would be some of the largest amounts of 
such dust observed in SNe to date, providing significant 
credi bility to SN dust models (e.g.. iNozawa et al]|2003L 
120081 ). Before any conclusions are drawn, however, we 
must first disentangle the composition, origin, and heat- 
ing mechanism of the dust. 

3. ANALYSIS: DUST ORIGIN AND HEATING MECHANISM 

The dust can have several different origins and heating 
mechanisms. It may be newly formed or it may have been 
in place at the time of the explosion. If newly formed, 
the d ust may have condense d in the expanding SN ejecta 
(e.g., Elmhamdi et al. 2004) or in the cool, dense shell of 



post-sh ocked gas lying be t ween the forward a nd reverse 
shocks (jPozzo et al.l l2004t iSmith et al.ll2008bD . In both 
cases, several heating mechanisms are possible, including 
radioactivity, radiative heating by the optical emission 
from circumstellar interaction, and collisional heating by 
hot gas in the reverse shock. 

Alternatively, pre-existing dust may be collisionally 
heated by hot, shocked gas or radiatively heated by ei- 
ther the peak SN luminosity or the late-time optical 
emission from circumstellar interaction. In the radia- 
tive heating case, the dust reprocesses the optical light 
and thermal radiation persists until the dust grains cool 
sufficiently. If this pre-existing dust is distributed in a 
shell with a light-crossing time greater than the duration 
of the optical emission, an "IR echo" is evident due to 
light-travcl-time effects. Multiple scenarios can also con- 
tribute to the late-time IR flux, as in the cases of SN e 
2004et (IKotak et all [2009]). 200 4di (jMeikle et al.l[201ll) . 
and 2006ic (|Mattila et al.ll2008h . 

With a detailed analysis of SN 2005ip. lFox^t~aII ([2TJl~0f ) 
illustrate how to use late-time mid-IR observations, along 
with optical photometry and spectra, to disentangle the 
various dust models. First, a few variables must be de- 
fined. The blackbody radius 



rbb = 



( L 



bi.) 



V4™r b 4 b 



1/2 



(2) 



where a is the Stefan-Boltzmann constant. The calcula- 
tion of the blackbody radius assumes an optically thick 
dust shell. While our calculations assume the case of op- 
tically thin dust, the blackbody radius sets the minimum 
shell size of an observed dust component. 
The shock radius 



v s t, 



(3) 



for a constant shock velocity v s , defines the maximum 
radius that the forward shock can travel in a time t. 
The shock velocity corresponds to the optical spectral 
line widths. Three width components are typical for 
SNe Iln (see Figures and El as well as 32. 2j) : narrow 
(< 700 km s^ 1 ), intermediate (~ 1000-5000 km s" 1 ), 
and broad (~ 5000-15,000 km s _1 ). The narrow lines, 
which give the Type Iln subclass its name, originate in 
the slow, dense, pre-existing circumstellar environment 
when excited by X-ray and ultraviolet (UV) emission 
generated by the forward shock. The intermediate lines 
correspond to the decelerated forward-shock front as it 
passes through the dense clumps in the CSM. The broad 
lines, which set the shock velocity, arise from the unin- 
hibited forward shock arising from the rapidly expanding 
SN ejecta. Over time, the dense CSM decelerates the en- 
tire ejecta. Typical forward shock velocities in SNe Iln at 
late times are v s « 5000 - 6000 km s" 1 (e.g ., ISmith et al.1 
[2W7l[Trmidle et al.ll2009t iRest et al.1120111) . 

Table [7] lists the associated blackbody and shock radii 
for each SN at the epoch of the Spitzer observation, as- 
suming a shock velocity of v s = 5000 km s _1 . Notice 
that unlike the shock radius, the blackbody radius is not 
a function of the observation epoch to the extent that 
the luminosity is constant. 

Other useful parameters are the peak supernova UV- 
optical luminosity, L pe ak, and the late-time optical lu- 
minosity, L opt . Dominated by shock heating of the stel- 
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Figure 6. The temperature, mass, and IR luminosity of the observed late-time dust component of each SN, as listed in Table [6] for 
a = 0.1 fim grains of bo th compositi ons. Downward arrow s indicate upper limits. Also plotted is the near-IR evolution of the Type Iln SNc 
2005ip l|Fox et alJI200Sh and 1995N HGcrardv et al.ll200^ . The number of SNe that exhibit late-time IR emission drops dramatically after 
~ 1000 days, and the trend appears to closely follow the observed decline of SN 1995N. Generally, the dust masses are M d > 10" 3 M Q , 
and many are even > 10 -2 Mq. Ultimately, these variables can be used to disentangle the composition, origin, and heating mechanism 
of the dust (see Q. As discussed in [T4] many of these results are consistent with LBV progenitors. 



lar envelope, the peak optical luminosity can vaporize 
all pre-exist ing dust within the vaporization radius (e.g., 
lDwea[l985lh 



' cvap 



7 pcak 



i67tctT c 4 (g) 



1/2 



(4) 



where T ovap is the vaporization temperature of the dust, 
and (Q) is the Pla nck-averaged val ue of the dust emis- 
sivity. Figure 8 of iFox et alj (|2010f ) shows the relation- 
ship between the peak luminosity and vaporization ra- 
dius. Given the limited number of published observa- 
tions for each SN during the discovery period, we obtain 



our best estimates of the peak luminosities from the well- 
known SN archival website, www.astrosurf.com/snweb2. 
The maximum reported luminosities for each supernova 
are in the range 3 x 10 8 < i pC ak < 4 x 10 9 L Q . Table 
[7] lists the corresponding vaporization radii, all of which 



are within 0.0065 < 



cvap 



< 0.02 ly. 



Whether a dust shell is located at or beyond r cvap , 
light-travel-time effects cause the thermal radiation from 
the dust grains to reach the observer over an extended 
period of time, thereby forming an "IR echo" (e.g.. 
Bode fc Ev ans 1985 |Pwe5[l983l: lEmmering & Chevalier! 



1988T 



see M3.3I for further discussion). The IR echo dura- 
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Table 7 

Blackbody, Shock, Echo Plateau, and Vaporization 
Radii 



SN 


Epoch 




~ b 
' s 




r d 




(days) 


(iy) 


(iy) 


(iy) 


(iy) 


2005cp 


1523 


0.038 


0.066 


2.1 


0.012 


2005gn 


1479 


0.016 


0.064 


2.0 


0.013 


2005ip 


948 


0.031 


0.040 


1.3 


0.01 


2006jd 


1149 


0.067 


0.050 


1.6 


0.01 


2006qq 


1048 


0.048 


0.044 


1.4 


0.016 


2007rt 


780 


0.057 


0.034 


1.1 


0.02 


2008J 


593 


0.035 


0.026 


0.81 


0.02 


2008cg 


474 


0.036 


0.020 


0.65 


0.02 


2008cn 


386 


0.031 


0.016 


0.53 


0.011 


2008gm 


301 


0.004 


0.014 


0.41 


0.0065 


2008ip 


438 


0.008 


0.020 


0.60 


0.015 



a The blackbody radius given by Equation [2] 

b The shock radius given by Equation[3]for u s =5000 

km s — - 1 . 

c The IR echo radius given by Equation [5] 

d The vaporization radius is the radius at which 

the dust temperature equals the dust vaporization 

temperature. 



SN 2006jd 




km s 1 



Figure 7. Spectroscopic evolution of the Ha line in SN 2006jd at 
late times. The underlying continuum has been subtracted. These 
spectra show no evidence of absorption. In fact, they exhibit an 
excess likely due to the [N II] A6583 line. 



tion t ec h defines the echo radius 



(5) 



Table [7] also lists the associated echo radii for each SN, 
where t cc h is a lower limit determined by the epoch of 
the latest Spitzer observations. 

At later times, once the radioactive component fades, 
interaction between the reverse shock can generate X- 
rays and UV radiation (e.s;.. | Draine fc Woods! H99l 
lChugailH99l iChevalier fc Franssonlll994l ). In the dense 
CSM, this radiation may be reprocessed and emitted 
at optical wa velengths, L nn t. For e x ample, for both 



SNe 2005ip (llmmler fc Poolevl EHM iFox et alJ EfflOl 
ISmith et al.ll2009cn and 1995N (IGerardv et al.l^OOa Tthe 
observed late-time X-ray luminosity is lower than the op- 
tical/IR, and for SN 1995N, the UV accessible to the 
Hubble Space Telescope has a power comparable to the 
optical (jFransson et alJ f2002f) . While late-time optical 
observations do not exist for most of the SNe in this 
paper, the late-time optical observations of SNe 2007rt 
( Trundle et all [2001) . 2006jd (sec Table O, and 2005ip 
Smith et al.ll2009q) serve as useful references. 



3.1. Newly Formed Dust 

New dust may condense in either the slowly moving 
ejecta or in the cool, dense shell behind the forward shock 
(and in front of the reverse shock) . Since the dust will 
preferentially absorb emission from redshifted (far-side) 
material, the lack of attenuation in the red wing of the 
optical emission lines provides a straightforward method 
for testing whether the dust is newly formed. (Due to 
the SN geometry, dust in the slowly moving ejecta closest 
to the core absorbs redshifted emission from the broad, 
intermediate, and narrow lines, while dust in the cool 
dense shell beyond the fastest ejecta absorbs redshifted 
emission only from the intermediate and narrow compo- 
nents.) 

Typically, such emission-line evidence for dust forma- 
tion occurs at later epochs, once the ejecta expand and 
become optically thin. For the first 40-50 days, larger 




4-1 0" 5-10* 



Figure 8. Spectroscopic evolution of the broad component of the 
Ha line in SNe 2008cg and 2008 J. An approximate fit of the under- 
lying continuum has been subtracted. It shows increased suppres- 
sion of the red wing relative to the blue wing with time, especially 
in SN 2008cg, which offers evidence of new dust formation in the 
rapidly expanding ejecta. 



optic al depths in th e ejecta suppress the redshifted emis- 
sion (]Chugaill2001l) so that attenuation in the red wing 
is little indication of newly formed dust. ISmith et alJ 
(|2009d ) spectroscopically confirm new dust formation in 
SN 2005ip as early as day 60 post-discovery, similar to 
the "Type Iln-likc" SN 2006jc, in which new dust was 
confirmed as early a s 51 days after t he pe ak brightness 
(|Smith et al.ll2008bD . iTrundle et ail (|2009D also provide 
spectroscopic evidence of dust formation for SN 2007rt. 
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Day 10 , 
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Table 8 

Predicted Dust Masses From Shock Heating 
(v B = 15,000 km s _1 ) vs. Observed 



Figure 9. Spectroscopic evolution of the intermediate component 
of the Ho line in SNe 2008cg and 2008 J. An approximate fit of 
the underlying continuum has been subtracted. The figure reveals 
an increased suppression of the red side relative to the blue side 
with time. The simultaneous absorption in the broad component 
(Figure [Sj , however, suggests that the dust forms in the ejecta. 



Of the remaining SNe in this sample, late-time spectra 
exist only for SNe 2006jd, 2008cg, and 2008J (see Table 
[4|) . Figures and [9] plot the Ha line from these op- 
tical spectra with the underlying continuum subtracted, 
and the velocity scale chosen with the narrow Ha emis- 
sion feature at v = km s -1 . In some cases, the nega- 
tive (blueshifted) velocity axis is mirrored (dotted line) 
on the positive (redshifted) side to provide a qualitative 
comparison between the redshifted and blueshifted com- 
ponents. Figure [7j shows no evidence of absorption in 
SN 2006jd (in fact, it exhibits an excess likely due to the 
[N II] A6583 line). 

Alternatively, Figures [5] and |H] show the broad (~ 
10 4 km s -1 ) and intermediate (~ 10 3 km s _1 ) compo- 
nents of the Ha line of SNe 2008cg and 2008J at late 
times. Both SNe show some degree of increasing absorp- 
tion in the red wing relative to the blue wing from days 
>63 and 224, respectively. Simultaneous absorption in 
both the broad and intermediate components may be 
evidence for newly formed dust in the ejecta, but other 
origins exist. For example, since elec tron scatterin g may 
be responsible for line broadening ([Chugail |2001D . it is 
plausible that this scattering could lead to an asymmet- 
ric line profile. Alternatively, the scattering opacity of 
dust may produce asymmetric line profiles even when 
the dust is entirely outside the shock. Even if the dust 
were newly formed, however, the masses listed in Table 
M(M d < 10" 2 M ) are less than in SN 2005ip and much 
lower than th e M d = 0.1-1 pr e dicted by the SN dust 
models fe.g.. iNozawa et all 120031 I2008D . Furthermore, 







Graphites 


Silicates 


SN 


Predicted 


Observed 


Observed 




Md (M ) 


M d (M ) 


M d (M ) 


2005cp 


4.9e-3 


6.8e-3 


1.6e-2 


2005gn 


4.6e-3 


1.9e-3 


3.2e-3 


2005ip 


1.9e-3 


4.3e-2 


9.3e-2 


2006jd 


2.8e-3 


3.0e-2 


5.0e-2 


2006qq 


2.3e-3 


1.6e-2 


2.8e-2 


2007rt 


1.3e-3 


2.2e-2 


3.8e-2 


2008J 


7.4e-4 


8.7e-3 


1.5e-2 


2008cg 


4.7e-4 


1.3e-2 


3.0e-2 


2008en 


3.1e-4 


6.7e-3 


l.le-2 


2008gm 


1.9e-4 


6.5e-5 


1.4e-4 


2008ip 


4.0e-4 


4.2e-4 


7.3e-4' 



the relative contribution of this newly formed dust to 
the overall mid-IR flux is not immediately clear. Along 
with the fact that the remaining SNe in this sample do 
not have late-time spectra, we must therefore consider 
the pre-existing dust models for all of these SNe. 

3.2. Shock Heating 

In the shock heating scenario, hot electrons in the 
post-s hocked gas collisiqnally heat pre-existing dust 
grains (|Draine fc SalDeterll97aiDrainell98ltlDweyi987l: 



IDwek et al.l 120081) . IFox et all (|2010fl derive the esti 
mated mass of dust heated by the forward shock, which 
serves as a useful consistency check when compared to 
the observed dust mass. For the forward shock at time t, 
the mass of the shocked gas is 



M g (M ) 0.28 



15,000 km s" 



yrj V/ in V 



(G) 



for a hydrogen atomic mass toh, shock velocity i> s , and 
grain size a, provided that the sputtering time scale 
Tsputt < (i, Tcooi), where t coo i is the radiative cooling time 
scale. The dust mass is then given by assuming a Galac- 
tic dust-to-gas mass ratio, Z d = M d /M g « 0.01. 

Table [8] compares the expected dust masses for the 
fastest shocks (v s = 15,000 km s _1 ) to the observed 
masses listed in Table [S SNe 2008ip, 2008gm, 2005gn, 
and 2005cp have dust masses consistent with the mod- 
els, but these cases all assume a forward shock traveling 
at a constant speed of v s w 15,000 km s _1 throughout 
the entire expansion without experiencing any decelera- 
tion. This is not likely for the two events from 2005, but 
it may be possible for the more recent SNe 2008ip and 
2008gm. The predicted dust masses for the remaining 
SNe, however, are more than an order of magnitude lower 
than observed. Furthermore, the dust-to-gas mass ratio 
in the stellar winds is likely lower than that of the ISM 
([Williams et al.l 120061 : ITvezic fc ElitzTuI 120101 ) . A lower 
ratio would predict even less dust. These results likely 
rule out shock heating for most of these events. 

3.3. Radiative Heating 

Optical emission from the SN may radiatively heat a 
shell of pre-existin g dust at a radius r to a tempera- 
ture Td. I Fox et aT] (|2010l) outline three basic scenarios. 
For the first two scenarios, the peak S N luminosity h eats 
the dust and forms an IR echo (e.g., [Wright 1983)- In 
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the first scenario, the dust may be distributed spheri- 
cally symmetrically around the star at the time of the 
explosion (most likely formed by a steady wind from 
the progenitor). The SN peak luminosity then vapor- 
izes the dust out to the vaporization radius r cvap , and 
heats the inside of the remaining dust shell to a tem- 
perature roughly equal to the vaporization temperature 
(Td re T cvap ). Alternatively, the dust may be distributed 
in a shell at a radius r > r ovap (most likely formed by a 
progenitor eruption hundreds to thousands of years prior 
to the SN). The SN peak luminosity heats the inside of 
this shell to a temperature Td < T cvap , where the rela- 
tionship between the two temperatures is generally given 



Ta — T P , 



L 



opt 



(7) 



where L opt is the observed SN optical luminosity and 
-^cvap is the luminosity required to vap orize the dust a t 
the given radius. In many cases 

(e.g.. lFoxet~aI|[2009[) . 

the resulting IR luminosity echo plateaus on year-long 
time scales, corresponding to the light travel time across 
the inner edge of the dust shell, where the shell radius is 
given by Equation [5] 

In the final scenario, the dust may be distributed in 
a shell at radii between the vaporization and echo radii, 
T'cvap < t < ^cch (formed by either a steady wind or pro- 
genitor eruption). The late-time optical luminosity gen- 
erated by circumstellar interaction, L opt , continuously 
heats the dust shell. This scenario is not so much a 
traditional IR echo as it is a reprocessing of the optical 
emission by the dust 03 If the circumstellar interaction 
occurs on a time scale greater than the light travel time 
across the dust shell (which we assume it does in these 
cases), the shell radius does not set the IR echo plateau 
length. The observed flux therefore accounts for the en- 
tire shell. 

For a simple light-echo model dominated by flux from 
only the innermost, warmest dust at a single tempera- 
ture, Td, balancing the energy absorbed and emitted by 
the dust grains gives the relationship between the SN 
optical luminos ity and the dus t temperature and inner 
shell radius as (jFox et al.ll201Qf ) 



T 64 2 ^4 



SN 



/ B v {T d ) K {v)dv 
J B u (T s ^)Q &hs (v)dv' 



(8) 



for a dust bulk (volume) density p and an effective SN 
blackbody temperature Tsn- When using this equation, 
we generally assume a = 0.1 //m grains and a SN tem- 
perature Tsn ~ 10,000 K, although the results are fairly 
insensitive to the choice of temperature. The luminosity 
is treated as a central point source, assuming the emit- 
ting region is internal to a spherically symmetric dust 
shell. 

Using Equation Figure [TO] plots contours of the ex- 
pected dust temperature as a function of shell radius r 
for various optical luminosities. We choose luminosities 
10 7 < L pt < 10 11 Li©, which span the observed peak 



10 Some authors re fer to this as a circumstellar interaction echo 
IIGerardv et al. 2002), but the reader should note that the heating 
mechanism is radiation from the shocks, rather than direct colli- 
sional heating by the shocks themselves. 



and late-time optical luminosities for these SNe. Over- 
plotted are the observed dust temperatures and radii de- 
rived from the Spitzer data (see Tables [6] and [7]) . Only 
dust with graphite composition is considered, as silicate 
grains require much higher luminosities than observed 
(> 10 11 L Q ) to heat dust grains of comparable sizes to 
the observed temperatures. 

Considering the first scenario above, the SN peak lu- 
minosity vaporizes all dust out to the vaporization ra- 



dius, 



' cvap s 



which is defined by the radius at which the 



dust temperature equals the dust vaporization temper- 
ature (for graphite, T cvap re 2000 K). None of the ob- 
served dust temperatures in Table [5] even approach the 
vaporization temperature of graphite. Furthermore, we 
can rule out this scenario by comparing the expected 
vaporization radii to the observed duration of the IR 
echos. Table [7] shows that the vaporization radii are 
all r ovap < 0.02 ly, which is too small to produce an 
IR echo over the observed year-long time scales. Of 
course, the unobserved shock breakout in the minutes to 
hours follow ing the supernova can approach peak value s 
- 10 11 L© llSoderberg et a!Jl200l iNakar fc Saril l20"l0h . 
Figure 8b in lFox et alJ (|2010f > shows that such values can 
vaporize 0.1 /im dust grains out to r ovap re 0.2 ly. Still, 
these radii are too small to produce an IR echo over the 
observed year-l ong time scales. 

Figure 10(a) tests the second scenario by comparing 
the expected and observed temperatures assuming the 
dust lies at the echo radii r oc h defined by the IR echo 
plateau duration i cc h- At these radii, the SNe require 
large peak luminosities (> 10 10 L Q ) to power the ob- 
served echo. Not only are the required peak luminosi- 
ties significantly greater than the measured ones, but 
they are larger than thos e of the most lumin ous core- 
collapse events ever found (jQuimbv et al.|[2007t) . The un- 
observed shock breakout luminosities are indeed higher, 
but the breakout typically lasts only for about 30 min- 
utes, even in the case of an extended red supcrgiant 
progenitor. While the shock breakout luminosity may 
be large enough to heat the dust to the observed tem- 
perature, the total radiated energy from the breakout 
(-^breakout « 10 11 L x 30 min < 10 48 erg) is not 
nearly enough to power the observed echo (£7echo re 



10 85 L x 1000 d re 10 5U erg), assuming an IR lumi- 
nosity L\n re 10 8 5 L Q extending throughout the 1000 
day period (see Figure ■ Even if we take an echo last- 
ing for only 100 days, the total radiated IR energy is still 
too large to have been powered by the shock breakout. 
This fact makes unlikely any IR echo powered by the 
peak SN lumino sity from shock breakout. 

Figure 10(b) tests the final scenario, in which the op- 
tical luminosity generated by continuous interaction be- 
tween the forward shock and CSM offers an alternative 
heating mechanism. In this case, the dust is assumed 
to lie at the smallest radius allowed by the observations, 
defined by the blackbody radius, r DD - 

The observed dust temperatures for most of the SNe 
require optical luminosities 10 7 < L opt < 10 8 L Q . While 
such late-time optical luminosities are larger than ob- 
served in most SN subclasses, they are quite realistic 
for Type Iln events. For example, late-time observa- 
tions show L opt (S N 2005ip) = 10 7 8 L on days >900 
(|Smith et al.ll2009cfl . L opt (SN 2006jd) = 10 7 7 L on day 
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Figure 10. Testing the possibility of radiative heating of a pre-existing dust shell (only graphites are considered here). The dust 
temperature, Tj, is given as a function of radius, r, for contours of constant optical luminosity. Overplotted are the observed dust 
temperatures and radii. The arrows signify that these values are all upper limits given the lack of constraints in our fitting routines. In 
case (a), the dust shell lies at the echo radius r Bcll , defined by the observed light-curve duration, and is heated by the peak SN luminosity, 
which is typically in the range 10 9 < L pt < 10 10 Lq. The required luminosities to heat the dust are too high. In case (b), the dust lies at 
the minimum radius given by r bt> and is heated by optical emission generated by continuous circumstcllar interaction, typically observed 
to be 10 7 < L op t < 10 8 . Aside from a few caveats discussed in the text, the observed temperatures are consistent with the expected 
values given by the contours. 



1609 (Table Ell, and £ op t (SN 2007rt) « 10 8 L© on day 
562 (ITrundle et aLll2Tj09lr rr l . 

The dust temperatures for some of the SNe, particu- 
larly SNe 2007rt and 2006jd, do fall above the expected 
values for a luminosity L pt ~ 10 8 L Q . However, in these 
cases a second, hotter dust component may be contribut- 
ing to the mid-IR flux. Already, we know from the fits 
in Fi gureElthat SN 2006.jd has a hotter dust component. 
Also. ITrundle et all (|2009f ) show that new dust forms in 
SN 2007rt, although they do not measure the temper- 
ature. Similarly, Figures [5] and [§] exhibit evidence for 
new du st formation in SNe 2008cg and 2008J. As in SN 
2005ip (|Fox et al.H201Q[ ). this newly formed dust may be 
hotter than the mid-IR dust component. In these cases, 
the observed Spitzer temperatures would be only upper 
limits. The resulting lower dust temperatures and corre- 
sponding lower blackbody radii would place these points 
significantly closer to their luminosity contour lines. 

Late-time optical observations are not available for SNe 
2005cp, 2005gn, 2006qq, 2008J, 2008cg, 2008en, 2008gm, 
and 2008ip, but almost all of these points fall close to or 
along luminosity contours comparable to observed late- 
time optical luminosities of SNe 2005ip, 2006jd, and 
2007rt. In general, the data are fairly consistent with 
heating by optical emis sion generated by con tinuous cir- 
cumstellar interaction. iGerardv et alj (|2002[ ) reached a 
similar conclusion in their study of late-time emission 
from five SNe. Two of the Gerardy et al. targets, 
SNe 1999Z and 1997ab, are actually on our initial tar- 
get list. The fact that we failed to detect any late-time 
mid-IR emission comes as little surprise given that the 
forward shock would eventually destroy the pre-existing 
dust shells, a trend we observe in Figure [6] and discuss 

1 1 We reconcile the fact that our most recent optical observations 
do not detect SN 2007rt by noting that they were made nearly 500 
days following the original Spitzer observations. 



more below in []4] 

A simple calculation shows that the pre-existing dust 
shells in this model are optically thin at IR wavelengths, 
consistent with assumptions for deriving the dust mass 
in Equation [TJ For a dust shell given by a blackbody 
radius r^b, the optical depth can be written as 



Tbb 



J<:1 



, o '"avgi 



(9) 



where K avg = 435 cm 2 g 1 is the absorption coefficient 
for graphite av eraged over 1-15 /zm (see Figure 4 of 
IFox et al.|[20l0h . For each SN except SN 2005ip, r bb < 1. 
Using the a dditional s pectro scopic data for SN 2005ip 
provided bv lFox et al.l (|2010[ ). however, we derive an op- 
tical depth Tbb ~ 0.5. Furthermore, since the blackbody 
radii are only lower limits, these optical depths are all 
upper limits. 

4. TRENDS, PROGENITORS, AND OVERALL MODEL 

For the above scenario, in which UV and optical emis- 
sion from circumstellar interaction radiatively heats pre- 
existing dust, the observed dust mass accounts for the 
entire dust shell since the shell-radius light travel time 
(2r b b/c) is significantly smaller than the circumstellar 
interaction time scale given by the epoch of observa- 
tion (see Table [7]). Assuming a dust-to-gas mass ra- 
tio expected in the H-rich envelope of a massive star, 
Zd = Ma/Mg w 0.01, the pre-existing dust mass corre- 
sponds to a progenitor mass-loss rate 



M- 



Z d Ar 

Mm 



120 km s" 



0.05 ly 



(10) 

ArJ yr v ; 
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Table 9 

Mass-loss Associated with Radiatively Heated 
Graphite Dust Shell 



S.N 


M x (10Ar)/r t c 
(M S yr" 1 ) 


'uptior 

(yr) 


-^opt 

log (L / L ) 


2005cp 


2.2e-3 


30 


I .40 


2005gn 


7.5e-4 


25 


7.10 


2005ip 


1.8e-3 


235 


7.88 


2006jd 


2.8e-3 


106 


7.68 


2006qq 


9.2e-3 


17 


7.72 


2007rt 


1.2e-3 


181 


7.79 


2008J 


1.6e-3 


55 


7.43 


2008cg 


2.3e-3 


57 


7.59 


2008cn 


1.4e-3 


19 


7.38 


2008gm 


1.6e-4 


4 


6.43 


2008ip 


3.3e-4 


13 


6.75 


a The expected optical luminosity g 
forward shock, given by Equation 112 


enerated by the 



for a wind speed v w . The SNc in W2.2\ have narrow line 
profiles with FWHM in the range ~ 100-500 km s" 1 . 
Since these narrow lines originate in the slow, dense, 
pre-existing circumstellar environment, they are typi- 
cally assumed to directly correspond to the progenitor 
wind speed. Assuming a thin shell, Ar/r = 1/10, Table 
[9] lists the associated mass- loss rate for each SN. 

The properties in Table |H] are generally consis- 
tent with LBVs during their giant eruptions, which 
have wind speeds on the order of hundreds of 
km s - 1 and mass-loss rates M = 10~ 5 — 10~ 4 M© 
yr -1 ([Humphreys fe Davidson I1994D . but can reach 



much higher (jSmith fe~wockill2006l:ISmith et al.ll2lf51 . 
In fact, the total dust-shell masses are in the range 
10~ 5 < M d < 10- 1 M Q (see Table©, which, assum- 
ing a dust-to-gas mass ratio Z d « 0.01, are also con- 
sistent with the range o f LBV shell masses observed by 
iSmith fc Owockil (|2006f) . A LBV has already been di- 
rectly identified as the lik ely progenitor of SN 2005gl 
(|Gal-Yam fc Leonard! I2009D . and, indirectly, LBVs have 
emer ged as the likely progenitors of many other SNe Iln 
'e.g.. iChugai fc Danzigerf 119941 : iSalamanca et al 



Gal- Yam et all 120071 : ISmith et all 120071 l2008al 
Kiewe et al.ll2010l ). 



12002 



2009a: 



Table [7] highlights the fact that in almost every case 
(SNe 2008gm and 2008ip aside), the blackbody radius 
is larger than the vaporization radius by nearly a fac- 
tor of two (rbb > 2r, 



cvap J 



If the dust were associated 
with a steady wind, the resulting shell size would have 
been more comparable to the vaporization radius. In- 
stead, these results suggest that the progenitor may have 
undergone a period of increased mass-loss prior to the 
SN at a time ^eruption ~ fbb/v w , which corresponds to 
tens to hundreds of years (see Table [9]) . Such an erup- 
tive e vent is inferred in many other Type Iln events 



e.g.. IChugai fc Danzigen 119941: ISalamanca et~al] 12002 . 

" 20071: ISmith et all 120071 I2008al 12009a : 



Gal- Yam et all ._ 
Kiewe et aljl2010f v 



Alternatively, Figure 8b in lFox et all (|2010T l shows that 
an unobserved sho ck breakout peak lum i nosity of only 
~ 10 10 L Q (jSoderberg et all 120081 : iNakar fc Saril 
120 101) can vaporize 0.1 /an dust grains out to r eV a P ~ 
0.05 ly, which is much more within the range of the ob- 
served blackbody radii. In this case, there is little evi- 
dence for an eruptive event. Instead, the progenitor more 



likely underwent a period of extreme, but steady, mass- 
loss prior to the SN, thereby enshrouding the progenitor 
in a massive gas and dust shell. 

To test each possibility, we can compare the observed 
late-time optical luminosities to the expected values for 
various CSM densities. We use the expression (e.g., 
Chugai 1992) to write the expected shock luminosity as 
a function of progenitor wind density, w = M/v w , 



L 



opt 



w- 



(12) 



where v w is the preshock wind speed and e < 1 is the 
efficiency of converting shock kinetic energy into visual 
light. While the conversion efficiency varies greatly de- 
pending on shock speed and wind density, we assume a 
value e « 0.5, acknowledging this value may be high. The 
preshock wind speeds are given by the narrow line widths 



reported in ^2.21 (and v w ~ 300 km s for unreported 
lines) . 

Table [9] lists the expected optical luminosity produced 
by a v s 5000 km s _1 forward shock moving through the 
CSM assuming a mass-loss rate equal to the value derived 
above for each dust shell. Except for SNe 2008ip and 
2008gm, the predicted values are nearly equal to the ob- 
served late-time optical luminosities reported above for 
SNe 2005ip, 2006jd, and 2007rt. This consistency sug- 
gests the wind density throughout the CSM is roughly 
equal to the value calculated for the dust shell, thereby 
favoring the extreme and steady wind scenario above. 

The story starts to get interesting when we revisit Fig- 
ure Overall, a significant percentage of the SNe Iln 
from the past 10 years (^15%) show late-time mid-IR 
emission (and the percentage is over 50% when only con- 
sidering SNe with explosion dates within 5 years of the 
observations.) After about 1000-2000 days (3-5 years) 
post-discovery, the luminosity begins to steadily decline. 
This decline may occur because the forward shock ul- 
timately overruns and destroys the dust shell. Alter- 
natively, the radiative efficiency of the shock will drop 
when the cooling time becomes longer than the expansion 
time. In either case, the apparent consistency amongst 
the detected SNe suggests that in these instances, many 
of the progenitors may have experienced similar mass- 
loss histories, although additional epochs of data will be 
necessary to confirm this possibility. If indeed the pro- 
genitors are LBVs, which are associated with initial stel- 
lar masses M > 25 M Q (jSmith et all 120041 ) . this result 
could have significant implications on mass-loss proper- 
ties of the high mass, l ow-metallicity stars obs erved in 
the early universe (e.g.. ISmith fc Owockil l2006h . These 
conclusions are also particularly noteworthy for stellar 
evolution theorists, whose models do not p lace LBVs in 
the final pre-SN phase (TSchaller et al Jll992t lLangedll993l: 
Lani rer et al.lll99i IStothers fc Chip|ll996t Mae der et all 
20051: IMaeder fc Mevnetl 1200811 



Of course, LBVs are not necessarily the progenitors of 
all SNe Iln, and several caveats should be mentioned. 
Primarily, there is a significant diversity among the tar- 
gets in this sample. Not only can some of the SNe 
with late-time emission be explained with other mod- 
els (e.g., shock heating), but a large percentage of the 
sample does not e xhibit any late-time X -ray, optical, IR, 
or radio emission (|van Dvk et al.l[T996h . Of well-studied 
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SNe Iln, many show vary i ng degrees of circum stellar in- 
teraction (|Li et all l200a ITrundle eTaTl [20090. In the 
extreme case of SN 1986J, for example, late-time opti- 
cal emissi on is observed more tha n two decades post- 
detection (jMilisavljevic et alJl200ll . (The faded Ha line 
relative to the strong [O I] line, however, suggests that 
the bulk of the emission is generated by the reverse shock, 
not the forward shock.) Not only do the observations 
vary, but several other SNe Iln ha ve shown evidence for 
red supergiant (RS G) progenitors (jFransson et al.ll2002t 
iSmith et al.ll2009bh . Furthermore, other authors invoke 
several additional scenarios to explain the dense circum- 
stellar shells, including wind-blown bubbles, high-density 
clumps, ve ry short Wo l f-Ray et stages following the LBV 
outbursts (|Dwarkadasl 1201 If ), mass ejection caused by 
frictiona l heating of a comm on envelope in a binary 
system ijTutukov et alJ Il992f ). and even relic disks of 
proto-stellar m aterial left over from the star formation 
(|Metzgerl[20Tl . 

5. CONCLUSION 

The Spitzer survey in this paper targets 68 SNe Iln 
discovered in 1999-2008 without late-time mid-IR obser- 
vations, as well as the well-studied Type Iln SN 1997ab. 
Ten of the 69 SNe (~ 15%) display late-time mid-IR 
emission from warm dust. The late-time optical and 
infrared observations presented here constrain the var- 
ious origins and heating mechanisms of the warm dust, 
which arc summarized by Table 1101 Although some SNe 
show evidence of newly formed dust or shock heating of 
pre-existing dust, observations for almost every SN are 
consistent with the scenario in which pre-existing dust 
is radiatively heated by optical emission from late-time 
circumstellar interaction between the forward shock and 
dense CSM (see the third scenario of ^3 . 3|) . A similar 
conclusion was reached by iGerardv et al.l (|2002[ ) regard- 
ing several other Type Iln events that have since faded. 

Although the model has not yet been confirmed, a 
trend amongst the Type Iln subclass begins to emerge. 
The late-time IR emission associated with this subclass 
appears to originate from pre-existing dust expelled by 
the progenitor star, which ultimately forms a shell when 
the interior dust is vaporized by the SN shock breakout. 
With some basic assumptions, the dust-shell properties 
indicate that the progenitor winds and mass-loss rates 
are most consistent with the properties of LBVs. In al- 
most every case, the data point to extreme, but steady, 
mass-loss histories as the progenitors evolve along toward 
their stellar death. These results therefore contribute 
to the growing evidence that SNe are not likely major 
sources of dust, although their progenitors might be if 
the forward shock docs not ultimately destroy the dust. 

Still, several caveats exist that offer alternative expla- 
nations for the dense dust shells without invoking LBV 
progenitors. The advantage to the model in this paper 
is that it allows us make predictions. If the dust is dis- 
tributed in a shell, and the heating mechanism is indeed 
continuous radiation generated by the forward shock, we 
expect the shock to eventually overtake the shell radius 
and destroy the dust. As the dust is destroyed, the late- 
time mid-IR emission should begin to significantly de- 
crease. While we see this general trend in Figure [HI no 
single SN Iln has multi-epoch mid-IR data that reveal 
the SN "turning off." 



Table 10 

Evidence for Dust Origin and Heating Mechanism 





Newly 


Shock 


IR 


Shock 




Formed? 


Heating? 


Echo? 


Echo? 


2005cp 


no 


maybe 


no 


yes 


2005gn 


no 


maybe 


no 


maybe 


2005ip 


yes 


no 


no 


yes 


2006jd 


no 


no 


no 


maybe 


2006qq 


no 


no 


no 


yes 


2007rt 


yes a 


no 


no 


maybe 


2008J 


maybe 


maybe 


no 


maybe 


2008cg 


yes 


no 


no 


yes 


2008cn 


no 


maybe 


no 


yes 


2008gm 


no 


yes 


no 


yes 


2008ip 


no 


yes 


no 


yes 



a ITrundle eTaTl 120091) . 



We can calculate the approximate point in time when 
this will happen as f ~ rbi,/v s , where v s is the forward- 
shock velocity. For the observed blackbody radii, rbb < 
0.1 ly and a shock velocity v s « 5000 km s _1 , we would 
expect the shock to overtake the shell at a time t < 5 yr. 
Once the shock destroys the dust, the luminosity decline 
will occur on time scales consistent with the dust-shell 
radius (tdcciino = 2rbb/c ~ 3 yr), as the dust destruction 
on the near side of the shell will be observed before dust 
destruction on the far side. This time scale is consistent 
with the apparent decline (determined from the upper 
limits) in Figure [6] Of course, the blackbody radius as- 
sumes an optically thick shell and only sets a lower limit. 
The actual dust radii are likely much larger. In addi- 
tion, slower shock speeds or asymmetric shell geometries 
would result in varying time scales, but mid-IR monitor- 
ing of these SNe will help to constrain the various models 
and geometries. 
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